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ABSTRACT 

The debris disk around the Vega-type star HD 34700 is detected in dust ther¬ 
mal emission from the near infrared (IR) to millimeter (mm) and submm wave¬ 
length range. Also detected is a distinct set of emission features at 3.3, 6.2, 7.7, 
8.6, 11.3 and 12.7//m, which are commonly attributed to polycyclic aromatic 
hydrocarbon (PAH) molecules. We model the observed dust IR spectral energy 
distribution (SED) and PAH emission features of the HD 34700 disk in terms of 
porous dust and astronomical-PAHs. Porous dust together with a mixture of 
neutral and ionized PAHs closely explains the dust IR SED and PAH emission 
features observed in the HD 34700 disk. Due to the stellar radiation pressure and 
Poynting-Robertson drag together with the photodissociation of PAHs, substan¬ 
tial removal of dust and PAHs has occurred in the disk, and continuous replen¬ 
ishment of these materials is required to maintain their current abundances. This 
implies that these materials are not primitive but secondary products probably 
originating from mutual collisions among planetesimals, asteroids, and comets. 

Subject headings: circumstellar matter — infrared: stars — stars: individual: 
HD 34700 


1. Introduction 

Vega-type stars are main-sequence (MS) stars with infrared (IR) emission in excess of 
the stellar photospheric radiation (see Backman & Paresce 1993). This “excess” emission is 
attributed to the thermal emission of solid dust grains in a “debris” disk, an optically thin 
disk orbiting around the central star. In debris disks, the gas is tenuous, and dust grains 
are removed easily by radiation pressure (RP) and the Poynting-Roberson drag (Burns et 
al. 1979, Backman & Paresce 1993). To sustain a debris disk, dust grains need to be 
replenished continuously by collisions between planetesimals, asteroids, and comets. This 
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implies that grains in a debris disk are not the primordial materials left over from the 
protostellar molecular cloud but remnants of the planet formation process. 

HD 34700, also known as SAO 112630, is classified as a Vega-like star based on its 1R 
excess (e.g., see Walker & Wolstencroft 1988). It has a spectral type of GOIVe, an effective 
temperature of T e g « 6000K (see Schiitz et al. 2009 and references therein), and a total 
stellar luminosity of L* fa 20.4 L & (Acke & van den Ancker 2004). With a Hipparcos 
parallax of fa 0.86 ± 1.84 mas, its distance d is very uncertain, which makes the estimate of 
its age, mass, and size ambiguous. It has been suggested that d is at least ~ 100 pc (Torres 
2004), and its apparent proximity to Orion provides an upper limit of d ~ 430 pc. Assuming 
d fa 200 pc, Torres (2004) estimated the stellar mass to be M* fa 1.1-1.2 M e . 

The HD 34700 system is known as a multiple system consisting of a spectroscopic binary 
(Torres 2004) and two other faint stellar components which are at a separation of ~ 5'.'2 and 
~ 9 / .'2, respectively, for Component B and Component C (Sterzik et al. 2005). While the 
two faint components are classified as bonafide T Tauri stars (TTS) based on their strong 
Ha emission and Li I A6798 absorption, the nature of the spectroscopic binary is somewhat 
vague. Both the Ha emission and the Li I A6798 absorption (with an equivalent width of 
~0.6 and ~ 0.17 A, respectively) detected in HD 34700 are relatively weak (Arellano Ferro 
& Giridhar 2003, Sterzik et al. 2005), indicating that the age of HD 34700 cannot be a few 
mega-years but a few tens of mega-years (Torres 2004). Moreover, as we will show later 
in Section 0 that the lifetime of the dust grains in the disk around HD 34700 is shorter 
than the stellar age, therefore it is plausible that HD 34700 has a debris disk rather than a 
protoplanetary disk. 

The HD 34700 debris disk has been extensively studied observationally, including space- 
borne and ground-based broadband photometry of the IR “excess” emission from the near-IR 
to submillimeter wavelengths. Of particular interest is the detection of the prominent 3.3, 6.2, 
7.7, 8.6, and 11.3 ^m emission bands in HD 34700. These emission features are collectively 
known as the “unidentified infrared” (UIR) emission bands (see Tielens 2008, Peeters 2014, 
Joblin 2015) and are now generally attributed to polycyclic aromatic hydrocarbon (PAH) 
molecules (Leger & Puget 1984; Allamandola et al. 1985). They are ubiquitously seen in 
a wide variety of Galactic and extragalactic regions, including protoplanetary disks around 
Herbig Ae/Be stars and TTS and disks in transition from a primordial to debris phase (e.g., 
see Brooke et al. 1993, Siebenmorgen et al. 2000, Acke & van den Ancker 2004, van Boekel 
et al. 2004, Geers et al. 2005, Sloan et al. 2005, Weinberger & Becklin 2005, Furlan et al. 
2006, Habart et al. 2006, Keller et al. 2008). However, the detection of the PAH emission 
features in debris disks is relatively rare. To the best of our knowledge, they are reportedly 
seen in only a few debris disks (Coulson & Walther 1995, Sylvester et al. 1997, Smith et al. 
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2004). For example, Acke et al. (2010) found that 37 out of 53 Herbig Ae/Be stars show 
PAH features in their 1R spectra obtained with the Infrared Spectrograph (IRS) on board 
the Spitzer Space Telescope. In contrast, Chen et ah (2006) reported no detection of PAH 
emission in the Spitzer/ IRS spectra of 59 MS stars classified as Vega-like objects based on 
the IRAS 60 pm excess. 

PAHs play an important role in gas-rich planet-forming disks. They dominate the heat¬ 
ing of the gas in their surface layers. They are also diagnostic of the presence of very small 
grains in their surface layers, which in return are diagnostic of the settling and coagulation 
processes occurring in the disks that ultimately lead to the formation of planetesimals (see 
J. Y. Seok & A. Li 2015, in preparation). For gas-depleted debris disks, PAHs could poten¬ 
tially be used to probe their disk structures that are widely considered as signposts for planet 
formation. Due to their stochastic heating nature (see Draine & Li 2001), the PAH emission 
features are expected to be spatially more extended than the emission of large grains that 
are in thermal equilibrium with the stellar radiation. Therefore, debris disks could be more 
easily resolved at the PAH emission bands. 

The 3.3 /irn PAH C-H stretching feature was clearly detected in the ~ 1.9-4.2 /irn spec¬ 
trum of HD 34700 obtained with the medium-resolution spectrograph SpeX at the NASA 
Infrared Telescope Facility (IRTF; Smith et al. 2004). The 6.2 and 7.7pm C-C stretching 
features as well as the 11.3 pm out-of-plane C-H bending feature were detected in the ~ 2.5- 
12 pm spectrum of HD 34700 obtained by the spectrophotometer ISOPHOT on board the 
Infrared Space Observatory (ISO) albeit its limited spectral resolution of A/AA ~ 90 (Acke 
& van den Ancker 2004), and later by Spitzer/ IRS with higher sensitivities (Berne et al. 
2009, Schiitz et al. 2009). By combining the ibS'(3/TIMMIS2 ~8-13 pm N-band spectrum 
with the Spitzer/ IRS spectrum, Schiitz et al. (2009) also reported marginal detections of the 
8.6 pm C-H in-plane bending feature and the 12.7 pm out-of-plane C-H bending feature. 

With a distinct set of prominent PAH emission features detected, HD 34700 offers a rare 
opportunity to explore the nature and origin of PAHs in debris disks. In this work we model 
the dust and PAH emission of the HD 34700 debris disk, with an aim of exploring its dust 
and PAH properties. This paper is organized as follows. We describe the observational data 
of dust and PAH emission in Section [21 We present the PAH and dust models in Section [3J 
The results are summarized in Section [4] and discussed in Section [5] with special attention 
paid to the photophysics of PAHs. The major conclusions are summarized in Section [6] 
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Table 1. Photometric Data for the HD 34700 System 


Telescope 

Wavelength 

Flux Density 

Uncertainty 

/Mission 

w 

(Jy) 

(Jy) 


0.36 (U) 

0.239 

0.0088 


0.44 (B) 

0.578 

0.0267 

EXPORT 

0.53 (V) 

0.887 

0.0409 


0.69 (R) 

0.941 

0.0780 


0.83 (/) 

0.887 

0.0818 


1.235 (J) 

0.968 

0.0205 

2MASS 

1.662 (H) 

0.847 

0.0179 


2.159 (K s ) 

0.678 

0.0150 


3.4 

0.467 

0.014 

WISE 

4.6 

0.327 

0.005 


12 

0.422 

0.004 


22 

3.33 

0.02 


9 

0.51 

0.012 


18 

1.76 

0.03 

AKARI 

65 

12.5 

0.8 


90 

10.9 

0.6 


140 

5.34 

1.0 


12 

0.605 

0.0424 

IRAS 

25 

4.42 

0.309 


60 

14.1 

1.41 


100 

9.38 

0.938 


450 

0.218 

0.037 

SCUBA 

850 

0.041 

0.0024 


1100 

0.039 

0.013 


1350 

0.0117 

0.0011 


Note. — The EXPORT UBVRI, IRAS, and SCUBA fluxes are taken from Mendigutia 
et al. (2012). The 2MASS, WISE, and AKARI photometric fluxes are extracted from the 
online catalogs. 
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2. Data 


We compile all the photometric data for HD 34700 available in the literature. The 
UBVRI, IRAS, and SCUBA fluxes are taken from Mendigutia et ah (2012) and references 
therein. The near- to far-lR fluxes are extracted from the 2MASS All-sky Point Source Cat¬ 
alog^ the WISE All-sky Source Catalog^ the AKARI Infrared Camera (IRC) All-sky Survey 
Point Source Catalogue (Version 1.0, Ishihara et al. 2010), and the AKARI Far Infrared 
Surveyor (FIS) All-sky Survey Bright Source Catalogue (Version 1.0) E All the photometric 
fluxes and uncertainties reported in the literature are summarized in Table [1] For the spec¬ 
troscopic data of HD 34700, we adopt the ~ 1.9-4.2 /mi IRTF /SpeX spectrum of Smith et 
al. (2004) and the ~ 5-14/mi Spitzer /IRS spectrum of Schiitz et al. (2009). A combination 
of both photometric and spectroscopic data enables us to perform a comprehensive modeling 
of the observed spectral energy distribution (SED) to examine all the major PAH features 
from ~ 3 to ~ 14 /mi as well as the dust properties simultaneously. 


Figure[T](a) shows the broadband photometry of HD 34700 together with the IRTF /SpeX 
and Spitzer /IRS spectra. In addition, we plot the Kurucz stellar atmospheric model spec¬ 
trum for T e fj = 6000 K and a surface gravity of Igg = 4.0 (Kurucz 1979). To correct for 
the interstellar extinction, the knowledge of the visual extinction (Ay) toward HD 34700 is 
required. In the literature, a range of A v values has been reported (e.g., A v ze 0-0.68 mag, 
Coulson et al. 1998, van den Ancker 1998, Acke & van den Ancker 2004), but it is as uncertain 
as its distance. In modeling the observed SED, we set the distance d and visual extinction 
Ay as free parameters. We search for the best values that reproduce the photometric fluxes 
with the Kurucz model. It is found that d ~ 260 pc and Ay 0.5 mag closely reproduce 
the photometric fluxes in the ultraviolet/optical and near-IR bands. Therefore, we adopt 
d = 260 pc with Ay = 0.5 mag for further analysis in this workQ 


Figures [l](b),(c) focus on the 3.3 /mi PAH feature and the PAH features at longer wave¬ 
lengths, respectively. It is interesting to note that the 3.3 /mi aromatic C-H stretch feature of 
HD 34700 is accompanied by a weaker feature at 3.43 /mi (see Figure [0(b)), which is generally 
ascribed to aliphatic C-H stretch (Joblin et al. 1996, Li & Draine 2012, Yang et al. 2013). 


1 http://www.ipac.caltech.edu/2mass/releases/allsky 

2 http: / / wise2. ipac.caltech.edu/docs/releases/allsky 

3 http: //www.ir.isas.jaxa.jp / AKARI / Observation/PSC/Public 

4 If Ay = 0, the same Kurucz model with d ss 320 pc can also closely reproduce the photometric fluxes up 
to the J, H, K s bands but is deficient in the WISE 3.4 and 4.6 gm bands. If this true, the excess emission 
at the WISE 3.4 and 4.6 gm bands may indicate the existence of an inner, warm “zodiacal cloud” in the 
HD 34700 disk. 
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In addition to the 3.43 /irn C-H stretch, aliphatic hydrocarbon materials also have two C-H 
deformation bands (see Li & Draine 2012). The Spitzer /IRS spectrum of HD 34700 reveals 
the presence of these two features at ~6.89 and ~ 7.23/mi (see Figure [0(c)). Compared 
with the PAH features seen in the diffuse ISM (e.g., see Ticlens 2008), the 7.7 /irn feature is 
somewhat red-shifted and appreciably broadened, while the 8.6 /irn feature is much weaker. 
These characteristics belong to either class B or C following the PAH spectral classification 
of Peeters et al. (2002), which is typically found in planetary nebulae and isolated Herbig 
Ae/Be stars. 


3. Model 

To simultaneously model the dust emission and the IRTF /SpeX and the Spitzer /IRS 
spectra, we adopt the porous dust model of Li & Lunine (2003a) and the astro-PAH moel 
of Li & Draine (2001) and Draine & Li (2007). 


3.1. PAHs 


Following Li & Draine (2001) and Li & Lunine (2003a), we assume a log-normal size 
distribution for the PAHs in the HD 34700 disk. The log-normal size distribution function 
is characterized by two parameters, o,q and a , which determine the peak location and the 
width of the log-normal distribution, respectively: 


(in pah / d(i — 


a/tt/S a {1 - erf [in (aPf n H /a 0 ) /V2 a] } a 


exp < —- 


ln(a/a 0 ) 


a 


r PAH 

lor a > a min 

( 1 ) 


We set the lower cutoff of the PAH size as = 3.5 A corresponding to ~ 20 carbon atoms 
in a PAH molecule (i.e., Nq ~ 20), which is the minimum survival size for PAHs in the 
diffuse ISM (see Li & Draine 2001). We treat a 0 and o as free parameters. The total PAH 
mass will be derived from fitting the observed PAH emission features. 


For the PAH absorption cross sections, we adopt those of Draine & Li (2007). Two 
types of PAHs are distinguished: neutral PAHs and charged PAHs. For charged PAHs, we 
do not distinguish either cations from anions or multiply charged PAHs from singly charged 
PAHs. 


The astro-PAH model of Li & Draine (2001) and Draine & Li (2007) represents the 
absorption cross section of each PAH feature by a Drude profile, which is expected for 
classical damped harmonic oscillators (see Li 2009a). Drude profiles are characterized by 
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Fig. 1.— (a) Spectral energy distribution of HD 34700. The photometric fluxes are shown with symbols 
(see Table [l]). The near-IR IRTF /SpeX spectrum (Smith et al. 2004) and rnid-IR Spitzer /IRS spectrum 
(Schutz et al. 2009) are shown as the orange and green lines, respectively. The photometric uncertainties 
are overlaid although most of them are smaller than the size of the symbols. Our best-fit model spectrum is 
presented as the red solid line, which is composed of the stellar (cyan dashed line), PAHs (blue dotted line), 
and dust (magenta dashed-dotted line) components, (b)-(c) Same as panel (a), but enlarged views of the 
IRTF /SpeX spectrum and the Spitzer /IRS spectrum of HD 34700, respectively, with measured uncertainties 
overlaid (gray). 










































their peak wavelengths, widths, and strengths (i.e., the wavelength-integrated absorption 
cross sections). As mentioned earlier (see Section [2]), the PAH feature profiles detected in 
HD 34700 somewhat differ from those typically seen in the diffuse ISM. In modeling the PAH 
emission features of HD 34700, we slightly adjust the peak wavelength and width for several 
features, but with the strength fixed. We note that it is not unphysical to adjust the peak 
wavelength and width as long as the strength is unaltered. 

We will only model the major PAH features such as those at 3.3, 6.2, 7.7, 8.6 and 
11.3 /jrn. Nevertheless, the intensities of the minor, aliphatic C-H features at 3.43, 6.89, and 
7.23 /irn are measured and will be briefly discussed in Section 15.11 The aliphatic C-H features 
of HD 34700 and other disks will be fully investigated in a separate paper (J. Y. Seok & A. Li 
2015, in preparation). 

The PAH features shown in Figure |T] manifest that both neutral and ionized PAHs 
exist in the HD 34700 disk because ionized PAHs emit strongly at 6.2, 7.7, and 8.6/rm, in 
contrast, neutral PAHs emit strongly at 3.3 and 11.3/rm. The charge distribution of PAHs 
is determined by the balance between the photoionization and the electron recombination 
(Bakes & Tielens 1994, Weingartner & Draine 2001), which is proportional to U -y/T gas /n e , 
where U is the ultraviolet (UV) starlight intensity, T gas is the gas temperature, and n e is the 
electron density. Moreover, the PAH size also affects its charging. Following the approaches 
described in Li & Lunine (2003a), we calculate the photoionization rate (/q on ) and the electron 
recombination rate ( k iec ) for PAHs of a given size at a given distance (r) from the central 
star (see Appendices A and B in Li & Lunine 2003a). 

To calculate k rec , we need to know the electron density n e as a function of r, the distance 
from the central star. We assume the cosmic-ray ionization of H 2 to be the major source 
of the electrons in the HD 34700 disk. The electron density is derived from n e ps r*<r cr^h 2) 
where <^cr ~ 3 x 10” 17 s _1 is the cosmic ionization rate and r* ~ lOMyr is the stellar 
age. We obtain wh 2 from M(H 2 ), the total mass of H 2 in the HD 34700 disk. We adopt 
M(H 2 ) = 1502 M® after scaling the original value of M(H 2 ) « 80 M® at d ~ 60 pc of 
Zuckerman et al. (1995) to d ~ 260 pc, the distance adopted in this work. Assuming a 
vertical hydrostatic equilibrium for the disk (see Equations (B2)-(B4) in Li & Lunine 2003a), 
we derive % 2 (r), the spatial distribution of H 2 , and then n e (see Figure Ef a)). Also, we take 
T gas (r) « (i?*/2r) 1 ^ 2 T e g « 591 (r/AU)~ 1//2 K, where if* m 0.0194 AU is the stellar radius (see 
Figure |2)(b)). In Figure |2](d), we show the ionization parameter U a/ T gas /n e as a function 
of the distance r from the central star, where the starlight intensity, U, is in unit of the 
912 A-l fim local interstellar radiation held of Mathis et al. (1983; MMP83). 

In Figure |3](a), the photoionization timescales (ri on = l/Aw) and the electron recombi¬ 
nation timescales (r rec = l/k rec ) as a function of PAH size are depicted for given distances 
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Fig. 2. — Radial profiles of (a) the electron density, n e (cm -3 ), (b) the gas temperature, T gas (K), (c) the 
starlight intensity, U, in unit of the 912 A-l /im MMP83 ISRF (e.g., see Equation (1) of Li & Lunine 2003a), 
and (d) the ionization parameter U\jT Eas /n e as a function of distance, r (AU), from the central star in the 
HD 34700 disk. Since the ionization balance of PAHs is controlled by U\/T gas /n e (Bakes & Tielens 1994, 
Weingartner & Draine 2001), the trend seen in panel (d) depicts the overall charging of PAHs depending on 
the distance to the central star. 
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Fig. 3. — Upper panel (a) Photoionization (non) and electron recombination (r rec ) timescales for PAHs 
in the size range of 3.5 A < a < 25 A at r = 10 AU (dotted lines), 15 AU (dashed lines), and 25 AU (solid 
lines) from the central star. In calculating r rec , we estimate n e from M( H 2 ), the mass of H 2 determined by 
Zuckerman et al. (1995), assuming that the cosmic-ray ionization of H 2 is the dominant contributor to the 
electrons in the HD 34700 disk. Lower panel (b) ionization fraction (^>i on ) of PAHs as a function of size at 
r =10 and 25 AU. 
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from the central star (r = 10, 15, and 25 AU). Note that r = 25 AU is considered to be the 
cutoff distance of the PAH distribution in the HD 34700 disk since at r > 25 AU the dust 
temperatures drop to T < 120 K (see Section 0|) , so the dust willx acquire an ice mantle 
and PAHs will condense onto the ice mantle. As discussed in Li & Lunine (2003a), r ion 
decreases for larger PAHs since the ionization threshold is inversely proportional to the PAH 
size. Also, Tion decreases as PAHs get closer to the central star. This is because the starlight 
intensity is enhanced (i.e., U oc r -2 , see Figure 121(c)). 

In comparison with Ti on , r rec is found to be much shorter than r ion at r = 25 AU whereas 
T rec becomes comparable to Ti on at r ~ 15 AU for large PAHs. At r = 10 AU, PAHs with 
a > 6 A have r ion even shorter than r rec - This indicates that PAHs, in particular large PAHs, 
can attain an ionized state in the inner region while PAHs in the outer region are more likely 
to be neutral (or negatively charged). 

Using the photoelectric emission and electron capture rates of Weingartner & Draine 
(2001), we investigate the PAH charging process in the HD 34700 disk. At a given distance 
from the star, we calculate the ionization fraction (0 ion ) of PAHs as a function of PAH size, 
which is the probability of finding a PAH molecule in a nonzero charge state. As shown in 
Figure |3](b) for 0 i on (a) at r = 10 and 25 AU, small PAHs are more likely to be neutral, and 
large PAHs are more likely to be ionized. This is consistent with the comparison between 
the photoionization and electron recombination timescales (see Figure [3]( a)). Overall, 0i on 
decreases as the distance r increases, consistent with the variation trend of U^T gas /n e with 
r (see Figure 12(d)) which determines the charging of PAHs as a function of r. Roughly 
speaking, PAHs near the central star are fully ionized since r ion r rec even for the smallest 
PAHs. At ~ 15-25 AU, however, because of r ion > r rec for PAHs of a < 20 A (see Figure |3](a)), 
a fraction of them will be negatively charged, so 0i on will be nonzero even for the smallest 
PAHs at r = 25 AU (see Figure [3](b)). Therefore, we expect a mixture of neutral and charged 
PAHs (cations or anions) in the HD 34700 disk. Note that, following Li & Draine (2001), we 
do not distinguish cations and anions when we model the IR emission of PAHs. 


3.2. Dust 

The porous grain model of Li & Lunine (2003a) for protoplanetary and debris disks 
assumes that the dust consists of porous aggregates of amorphous silicate and carbonaceous 
materials, and the dust in regions cooler than ~ 120 K is further coated with H 2 0-dominated 
ices. We take the dielectric functions of “astronomical” silicate of Draine & Lee (1984) for 
the silicate component, amorphous carbon of Rouleau & Martin (1991) for the carbonaceous 
component, and “dirty” ice of Li & Greenberg (1998) for the ice component. Assuming 
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a spherical shape, the absorption cross sections of porous dust are calculated from Mie 
theory combined with the Maxwell -Garnett and Bruggeman effective medium theories (see 
Li 2009a). We adopt the mass mixing ratios of Li & Lunine (2003b) derived from cosmic 
abundance considerations. We first employ the Maxwell-Garnett effective medium theory 
to calculate the average dielectric functions for the ice-coated silicate subgrains and the ice- 
coated carbonaceous subgrains (see Equations (7), (8) of Li & Lunine 2003a) and then use 
the Bruggman effective medium theory to calculate the average dielectric functions for the 
porous dust aggregate (Equation (9) of Li & Lunine 2003a). The porosity (or fluffiness) of a 
grain is quantified by the fractional volume of vacuum (P). Following Li & Lunine (2003a, 
2003b), we adopt P = 0.90 for porous dust and P' ~ 0.73 for the ice-coated porous dust, 
and the choice of P will be discussed in Section [4] 


For the spherical dust, we adopt a power-law size distribution: dn/da cx a - ", where a 
is the spherical radius of the dust, and a is the power-law index. We take the lower cutoff 
size to be a m i n = 1 /jrn and the upper cutoff size to be a max = 1 cm. Li & Lunine (2003a) 
have shown that changing a m i n and a max does not significantly alter the model-fitting. This 
will be discussed in Section [H We vary a from 0.1 to 4.0 to search for the best fit. 

For the spatial (radial) distribution of the dust in the HD 34700 disk, instead of using a 
simple power-law with sharp inner and outer boundary cutoffs, we propose a more physical 
formula for the spatial distribution, which is defined as 


dn 

dr 




( 2 ) 


where r min is the inner boundary of the disk, which is set to be the distance where the dust 
temperature reaches ~ 1500 K and dust sublimates. This functional form has the advantage 
that on one hand, it resembles a power-law dn/dr cx r -7 at larger distances (r ;§> r m i n ), and 
on the other hand it peaks at r p = r min (/3 + 7) /q, unlike the simple power-law which peaks at 
r min . The latter is unphysical as it assumes that the dust piles up at r min where dust actually 
sublimates. We take r min = 0.3 AU since the equilibrium temperature of micrometer-sized 
dust reaches ~ 1500 K at r ~ 0.3 AU from the central star. We take the outer disk boundary 
to be r max = 1000 AU. Our model is not sensitive to r max , but r max = 1000 AU is a reasonable 
assumption in view of the distance to the nearest companion star in the system of HD 347000 
In summary, we have two free parameters, r p and 7, for the dust spatial distribution. 


5 The angular distance from the central star to the stellar component B is ~ 5'.'2 (Sterzik et al. 2005), 
which corresponds to ~ 1350 AU at a distance of d = 260 pc. 
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Fig. 4.— Log-normal size distribution for the PAH population. The best-fit model with ao = 3.5 A and 
a = 0.3 is compared with that of the interstellar PAH model (ao ~ 3.5 A and cr ~ 0.4, Li & Draine 2001). 
The size distribution is expressed by multiplying a 4 to show the mass distribution per logarithmic PAH 
radius. 


4. Results 

Figure [Tj shows our best-fit model that reasonably reproduces all the observational data. 
The best-fit model includes emission from a mixture of neutral/ionized PAHs and porous dust 
with a power-law size distribution index of a ~ 3.5 and a total dust mass of Md us t ~ 21.9 M®. 
The disk geometry is characterized with r p 100 AU and 7^2 (see Equation (T2T) ). The 
PAH component is determined by ao ~ 3.5 A and a ~ 0.3, and the total PAH mass is 
A^pah ~ 4.53 x 10 -7 (r/AU ) 2 M® if all PAHs are at a distance of r from the central star. We 
do not have any prior information how the PAHs are spatially distributed in the HD 34700 
disk. Due to their stochastic heating nature (Draine & Li 2001), the PAH emission spectral 
profiles do not vary with the starlight intensity. Therefore, the PAH emission features should 
essentially remain identical for PAHs at different distances from the central star, except that 
the flux level scales with r~ 2 . Consequently, the PAH mass required to account for the 
observed PAH emission features scales with r 2 . 

The PAH size distribution with the best fit parameters is shown in Figure [4] together 
with that of the interstellar mixture (i.e., ao ~ 3.5 A and a « 0.4, Li & Draine 2001 ). 
Compared with the diffuse ISM, the PAH emission features observed in the HD 34700 disk 
are in favor of smaller PAHs. 
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Fig. 5. — (a) Temperatures of porous dust (of porosity P = 0.90 for T > f20K and P' = 0.73 for 
T < 120 K) as a function of grain size (a) and radial distance (r) from the central star, (b) Temperatures of 
porous dust as a function of grain size at r = 20 AU (dotted lines), r = 25 AU (solid lines), and r = 100 AU 
(dashed lines). A sudden temperature drop around a ~ 30 fim is caused by the ice sublimation at T « 120K. 
For comparison, the gas temperatures T gas « 591 (r/AU) -1 ^ 2 K at the same distances are overlaid. 
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Our model yields a total dust IR flux of F\d\ ps 1.31 x 10 9 erg cm 2 s 1 , of 
which the PAH component contributes r\j 10%p| This results in a total IR luminosity of 
Lm ~ 2-75 L q at d — 260 pc and a ratio of the IR luminosity to the stellar luminosity of 
Lir/L* ~ 0.13. Previous studies estimated Lir/L*~ 0.14 (Sheret et al. 2004) or ~0.15 
(Sylvester et al. 1996), in close agreement with our estimate. 

The best-fit model predicts that the peak radius of the dust spatial distribution in the 
HD 34700 disk is located at r p « 100 AU from the central star. The peak radius r p derived 
here can be compared with the inner radius (i.e., r min ) of the disk models proposed in the 
literature which assume a simple power-law dust spatial distribution and a peak distribution 
at f mm- Sylvester et al. (1996) and Dent et al. (2005) derived r min ~ 24-46 AU (at d ~ 55 pc) 
and < 30 AU (at d & 125pc), respectively. If we scale them to those at d — 260pc, our 
estimate is consistent with the former but is slightly larger than the latter. 

Using the power-law size distribution and the spatial distribution of the dust (see Equa¬ 
tion fl2ll). we estimate the total dust mass in the disk to be Md us t ~ 21.9M®. There have 
been various estimates of the dust mass of the HD 34700 disk (e.g., see Zuckerman et al. 1995, 
Coulson et al. 1998, Sheret et al. 2004, Dent et al. 2005). Although these previous studies 
adopted different distances to HD 34700, ranging from d « 60 pc to d > 180 pc, the derived 
dust mass ranges from Md us t ~ 7 to ~30M® after scaling to d = 260 pc (i.e., Md ust oc d 2 ), 
which overlaps with our estimate. 

We calculate the equilibrium temperature, T(r, a), of porous dust as a function of the 
radial distance r from the star and grain size a (see Figure (5](a)). The temperature decreases 
with increasing grain sizc@ and distance from the star. Figure [5](b) shows the equilibrium 
dust temperatures as a function of dust size at specific distances (r = 20, 25, and 100 AU). 
At r = 25 AU, we see a sudden temperature drop around a ~ 30 /mi. This is caused by 
the ice sublimation at T ^ 120 K: porous dust larger than a ~ 30 /mi attains T ss 120 K 
at r « 25 AU. When the silicate and carbonaceous subgrains acquire an ice mantle at 
T < 120 K, their aggregates will be less absorptive and therefore attain a lower equilibrium 
temperature. 

The location where dust grains have low enough temperatures to capture volatile molecules 
to form ice mantles is called the “snowline.” The snowline plays an important role in the 
PAH chemistry since free-flying PAHs will condense onto the ice mantles. As shown in Fig- 


6 In contrast, in the Galactic diffuse ISM the PAH component accounts for ~ 1/3 of the total IR luminosity. 

'For large grains with a > 100/rm, they emit like a blackbody and their temperatures are independent 
of grain size. 
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ure [5f a), the transition from non-icy grains to icy grains occurs at different distances from 
the star for grains of different sizes: for smaller grains, the “snowline” is further away from 
the star, and by r « 80 AU all grains become ice-coated. For simplicity, in the following, we 
adopt r — 25 AU as the snowline of the HD 34700 disk. 

Although the observed SED is closely reproduced by our best-fit model, there could exist 
some flexibilities for the model parameters. We briefly describe here the robustness of our 
model and refer to Li & Lunine (2003a), who explicitly discussed all the parameters of the 
porous dust model for disks (see their Section 4.1 for details). First, we examine the adopted 
dust parameters a m i n , a max , and P while other parameters remain unchanged. We differ a m ; n 
between 0.1 /irn and 10 /im. While models with a large a m i n (~ 10 /irn) emit too little in the 
mid-IR at A < 40 /mi, those with a small a min (~ 0.1 /mi) deviate the fit at the mid- and 
far-IR. Similarly, by varying a max from 100 /irn up to 1 cm, we find that sufficiently large 
grains (> 0.1cm) are necessary for reproducing the observed subnnn emission. The models 
with 1000 /i,m < a max < 1 cm produce almost identical SEDs. As a max decreases, the total 
dust mass (M<j ust ) decreases by a factor of ~ 3 at most. The porosity P is degenerate with 
a min and r min since decreasing a min or r min or increasing P has the same effect of increasing 
the dust temperature. However, a porosity in the range of 0.80 < P < 0.90 is expected for 
dust aggregates formed through coagulation as demonstrated both theoretically (Cameron & 
Schneck 1965) and experimentally (Blum & Wurm 2008). Moreover, a porosity of P ~ 0.90 
is consistent with the mean mass density of cometary nuclei for which the ice-coated dust 
aggregates are plausible building blocks (see Greenberg & Li 1999). It is worth noting that 
the fluxes calculated from models with 0.8 < P < 0.95 for the HD 34700 disk agree with 
each other within < 25%. 

We also examine the dust spatial distribution parameters r p and 7 which are not readily 
constrained due to the lack of spatial resolution of the disk. We find that models with small 
r p (< 80 AU) would emit too much around A ~ 10 /irn due to the silicate emission feature, 
while models with large r p (> 200 AU) fail to explain the far-IR/submm emission. The 
80 AU < r p < 150 AU models produce reasonable fits, with Md ust varying by a factor of ~ 2 
at most. We also End that varying 7 (1.0 < 7 < 3.5; combined with r p ) does not affect the 
model-fit significantly. In summary, although different combinations of model parameters 
may be able to explain the observed SED, the underlyling results do not change. The 
estimate of Md us t is most vulnerable (by a factor of ~ 3 at most), while the other estimates 
(e.g., Mrp and Mpr in Section [575]) vary within ~20%. 
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5. Discussion 

5.1. PAHs in the HD 34700 Disk 

As mentioned in Section EO due to the single-photon heating nature of PAHs, the total 
mass of PAHs in the disk (M P ah) required to reproduce the observed PAH emission is simply 
scaled by U, the starlight intensity. However, we cannot constrain U (which is inverse 
proportional to r 2 , i.e., U(r) oc (r/AU)~ 2 ) since the HD 34700 disk has not been spatially 
resolved and hence we do not have any prior knowledge of the spatial distribution of PAHs. 
One may derive M P ah ~ 2.83 x 10 ~ 4 M® by locating all PAHs at r = 25 AU0 the location 
of the snowline in the HD 34700 disk. Beyond the snowline, grains start to be ice-coated 
(see Figure [5]), and the free-flying PAHs will condense onto the ice mantles and will not be 
excited to emit at the observed PAH features. 

Even with Mpah ~ 2.83 x 10 _ 4 M® for r = 25 AU, the PAH-to-dust mass ratio is 
only MpAH/Md us t ~ 1.3 x 1O _5 0 much smaller than that of the Galactic diffuse ISM 
(MpAH/Md us t ~ 0.05, Li & Draine 2001, Draine & Li 2007). The deficit of PAHs rela¬ 
tive to dust has previously been reported for debris disks (e.g., see Li & Lunine 2003a, Thi 
et al. 2013). However, the derived PAH-to-dust mass ratio may not reflect the actual PAH- 
to-dust mass ratio in the HD 34700 disk as there could be a substantial amount of PAHs in 
the disk beyond the snowline. 

The minor features at 3.43, 6.89, and 7.23 /mi shown in the HD 34700 disk reveal that the 
PAH molecules very likely have aliphatic sidegroups, although the 3.43 /mi feature could also 
be due to anharmonicity and/or super hydrogenation of the aromatic C-H stretch (Barker et 
al. 1987, Bernstein et al. 1996). Let I^a and I 3 3 be the observed intensities of the 3.43/mi 
and 3.3 /mi emission features, respectively. For the HD 34700 disk, the IRTF /SpeX spectrum 
gives -I3.4/ 13.3 ~ 0.49. Let A 3 4 and A 3 3 be the band strengths of the aliphatic and aromatic 
C-H bonds, respectively. Yang et al. (2013) calculated A 3 . 4 /A 3.3 m 1.76. By assuming 
that the 3.43 /mi emission is exclusively due to aliphatic C-H (i.e., neglecting anharmonicity 
and superhydrogenation), we can place an upper limit on the ratio of the number of C 
atoms in aliphatic sidegroups to that in aromatic benzene rings: Y Cialiph /Ac )ar om ~ 0.3 x 


8 This means that most of the PAHs are concentrated within r ~ 25 AU corresponding to an angular size 
of 9 ps 0'/l at a distance of d ~ 260 pc. Future IR observations with a sub-arcsecond spatial resolution such 
as the Near Infrared Camera (NIRCam) onboard the James Webb Space Telescope will shed light on the 
spatial distribution of dust and PAHs in debris disks. 

9 In HD 34700, the fractional IR luminosity (Lpah/Lir ~ 10%, see Section 0]) emitted by PAHs is much 
higher than the PAH-to-dust mass ratio. This is because PAHs are closer to the central star and also much 
more absorptive than the bulk porous dust. 
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(I3.4/I3.3) x (A3.3/A3.4) ~ 0 .22. It is assumed here that one aliphatic C atom corresponds to 
2.5 aliphatic C-H bonds and one aromatic C atom corresponds to 0.75 aromatic C-H bond 
(see Li & Draine 2012). Similarly, we can place an upper limit of A^c, a ii P h/^c,arom ~ 0.14 
from the ratio of the observed intensity of the 6.89 //rn feature to that of the 7.7 /irn feature 
(J 6 . 89 / 77.7 ~ 0.064; see Section 3 of Li & Draine 2012 for details). 

It is not clear why PAHs are rarely seen in debris disks. This may be related to the disk 
structure and the stellar age. 

There exists evidence that most of the protoplanetary disks in which the PAH emission 
features have been detected have a flaring disk geometry (Meeus et al. 2001, Acke & van den 
Ancker 2004). A common interpretation for this is that flared disks intercept more stellar 
radiation than flat ones, especially at large distances from the central star; the PAH emission 
originates in the surface layers of a flared disk, where PAHs are directly exposed to starlight 
(Chiang & Goldreich 1997, Meeus et al. 2001, Acke & van den Ancker 2004, Habart, Natta, 
& Krugel 2004). It is not clear if this scenario applies to the HD 34700 disk as it is not 
spatially resolved. 

If the PAHs seen in debris disks originate from the outgassing of cometary bodies (i.e., 
the sublimation of the icy mantles coated on cometary grains onto which PAHs have con¬ 
densed during the protostellar dense cloud phase, see Li & Lunine 2003a), it is more likely 
to see PAHs in relatively young debris disks as they are subject to massive bombardment 
by cometary bodies (Chyba et al. 1994). With an age of a few tens of mega-years (Torres 
2004), the HD 34700 disk is probably still in the heavy bombardment phase as it lasted more 
than 500 mega-years in our early solar system (Chyba et al. 1994). 


5.2. PAH Destruction 

A PAH molecule is subjected to photodissociation upon absorbing an energetic photon. 
During this process, PAHs, preferentially small PAHs, might lose a hydrogen atom, a hydro¬ 
gen molecule, or an acetylene molecule (C 2 H 2 ), which eventually results in the destruction 
of a PAH molecule. Following Li & Lunine (2003a), we have calculated the photodestruction 
rates (fc des ) for the PAHs in the HD 34700 disk assuming that a PAH molecule is destroyed 
through the C 2 H 2 ejection. Figure E] shows the photodestruction timescales (r des = l/fcdes) 
at given distances (r = 10, 15, and 25 AU). During the lifetime of HD 34700 (r* ~ lOMyr), 
PAHs with a size of < 4.1 A (corresponding to PAHs with ~ 35 carbon atoms) are likely 
to be photodestroyed. Against the complete removal of such small PAHs in the HD 34700 
disk, these PAHs have to be replenished continuously, which is indeed required to account 
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Fig. 6.— Photodestruction timescales (rd es ) for PAHs at r = 10 AU (dotted line), 15 AU (dashed line), and 
25 AU (solid line) as a function of size. The age of HD 34700 (r+) is rather uncertain, ranging from ~ 3Myr 
up to several tens of Myrs. The age range (r* ~ 3 — 30 Myr) is designated by the shaded region with the 
representative value of r* ss 10 Myr plotted as a long-dashed line. 


for the observed PAH emission. Assuming the PAHs in the HD 34700 disk follow the same 
spatial distribution as porous dust (see Equation (121) ) but with a cutoff at r = 25 AU, we 
estimate the PAH mass replenishment rate to be Mpah ~ 4.5 x 10~ 9 M^yr^ 1 taking only 
photodestruction into account (see also Section 1531) . 


5.3. Radiation Pressure 

Small dust in a debris disk is continuously expelled outward due to the RP against the 
gravitational force (Backman & Paresce 1993). Following the approach described in Section 
4.3 of Li & Lunine (2003a), we derive the ratio of the radiative force to the gravitational 
force (/3rp) as a function of grain size.Both the RP and the gravitational force are inversely 
proportional to r 2 , as a result, /? RP is independent of the distance from the central star. 

We find that grains smaller than ~ 20-40 prn have /3 R p > 1 (see Figure [7]). This indicates 
that grains with a radius smaller than ~ 40 pm (~ 20 pm for icy grains) will be blown out due 
to RP. Moreover, assuming that the gas drag is not overwhelming which is true for debris 
disks, even larger grains (~ 40-80 pm) could be pushed outward (see Figure [TJ) because dust 
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Fig. 7. — Ratio of the radiative force to the gravitational force (/3rp, top panel) and the Poynting-Robertson 
drag timescale (tpr, bottom panel) for ice-free porous dust (P = 0.9: solid line) and ice-coated porous dust 
(P' ~ 0.73, dashed line) at r = 100 AU (note tpr oc r 2 ). The stellar age (r*) is marked in the same way as 
Figure [6] 
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particles released from parent bodies at periastron only need /3 rp = 0.5 for ejection (Burns 
et al. 1979). 

Unlike dust grains, the /3 rp ratio for PAHs is not sensitive to their sizes as they are in 
the Rayleigh limit (see Li 2009a). We estimate /3 rp ~ 35 for the PAHs in the HD 34700 disk. 


5.4. Poynting Roberson Drag 

Dust grains can migrate inward due to the Poynting-Roberson (PR) drag (Backman 
& Paresce 1993). We calculate the PR drag timescale (t pr ), which is a function of grain 
size and the distance (r) from the star. Note that, unlike the RP ejection, the PR drag is 
proportional to r 2 . 

At the peak distance of r p ss 100 AU, grains smaller than ~ 50-100 prn have a PR 
lifetime shorter than the age of HD 34700 (assuming r* ~ lOMyr; see Figure [3). The closer 
dust grains are located to the star, the bigger grains will be gradually accreted to the central 
star. For PAHs, we estimate tpr ~ 6000 yr for PAHs at r = 25 AU, which is independent of 
their sizes for the same reason as for /Irp. 


5.5. Dominant Processes for the Removal of Dust and PAHs in the HD 34700 

Disk 

Taking into account the RP ejection, the PR drag, as well as the photodissociation of 
PAHs, we examine the removal of dust and PAHs in the HD 34700 disk as a function of 
their sizes and distance from the star. The RP timescale in the HD 34700 disk is estimated 
to be t rp( r) !=s 912 (r/100 AU) 3//2 yr assuming that it is comparable to the local dynamical 
timescale 7d yn = 27r/f2(r), where Q(r) is the Keplerian frequency. For grains with /3 rp > 0.5 
(i.e., a < 80 pm), the RP timescale at a given distance is much smaller than the PR drag 
timescale (see Figure [7]). Therefore, for these grains the RP expulsion is the dominant removal 
process over the entire disk. For larger grains (i.e., a > 80 pm), if the PR timescale r PR (a, r) 
is shorter than the stellar age t*, the PR drag becomes an effective removal mechanism. At 
r = 100AU, this applies to grains of a < 100pm (see Figured). For grains at a larger 
distance from the star, the size range over which the PR drag is effective becomes narrower, 
and eventually, the PR drag no longer contributes to the dust removal beyond r ~ 150 AU. 

Since the overall dust removal in the HD 34700 disk is dominated by the RP ejection, 
we calculate the RP dust mass-loss rate for our best-fit model. Following Equation (20) in 
Li & Lunine (2003a), we derive Mrp ~ 2.05 x 10 -4 M e yr -1 . Hence, ~2000M® of dust 
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has been removed during the lifetime of HD 34700. For comparison, the PR dust mass-loss 
rate integrated over the dust size distribution and over the entire disk is estimated to be 
M PR »* 1.78 x 10 -7 M® yr -1 . 

In the same manner as for dust grains, we estimate the PAH mass loss rate due to the 
RP expulsion to be Mrp(PAH) « 2.92 x 10 _6 M®yr _1 . Note that PAHs are also subject 
to photodissociation (see Section 15.21) . We found that the RP timescale is much shorter 
than the photodestruction timescale except for the smallest PAHs (a ~ 3.5-3.6 A). The PAH 
mass loss rate via the RP expulsion (M R p) yields that ~ 30 M e of PAH molecules have been 
removed from the disk or have been replenished to maintain the current PAH abundance 
during the lOMyr lifetime of HD 34700. 

It is interesting to note that the PAH-to-dust mass-loss ratio, M RP (PAH)/M RP (dust) ~ 
1.4 x 10~ 2 , is about three order of magnitudes higher than the PAH-to-dust mass ratio of 
ATpah/A fdust ~ 1.3 x 10 -5 (see Section 15.11) . This is because only grains with a < 80/.an 
will be ejected by the RP expulsion (see Figure [7J) and with a best-fit dust size distribution 
of dn/da oc a~ 3 ' 5 (see Section QJ, the total dust mass in the HD 34700 disk is dominated 
by much larger grains. Thus, only a small mass fraction of the dust is subject to the RP 
expulsion, whereas PAHs of all sizes have been removed or replenished. 

The high removal rates of small dust and PAHs derived above require an efficient replen¬ 
ishment probably provided by the outgassing and collisions of planetesimals, asteroidal, and 
cometary bodies to contain a stable amount of small grains and PAHs in the disk (Backman 
& Paresce 1993, Wyatt 2008, Kiefer et al. 2014). In particular, as PAHs have been detected 
in interplanetary dust particles, comets, and meteorites (e.g., see Li 2009b for overview), 
cometary origin of PAHs in the HD 34700 disk is plausible. By systematically studying all 
debris disks with detected PAH features, we will (hopefully) be able to constrain the origin 
of PAHs and dust in debris disks (J. Y. Seok & A. Li, 2015, in preparation). 


6. Summary 

The HD 34700 disk is one of the few debris disks that show prominent PAH emission 
in their IR spectra. We have modeled the PAH emission of the HD 34700 disk as well as 
its entire SED from the near-IR up to submm/mm wavelength regime. The stellar RP and 
Poynting-Robertson drag quickly remove the dust and PAHs in the HD 34700 disk, implying 
that they must have been replenished continuously to maintain their current abundances. 

The major PAH features at 3.3, 6.2, 7.7, 8.6, and 11.3 //m (as well as the minor features at 
3.43, 6.89, and 7.23 /irn) are present in the IR spectra, indicating the presence of a mixture 
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of neutral and ionized PAHs (with aliphatic sidegroups) in the disk. By comparing the 
photoionization and electron recombination timescales, it is found that the PAHs in the 
inner disk are likely to be ionized while those in the outer disk are to be neutral or negatively 
charged. Also, at a given distance from the star, smaller PAHs are more likely to be neutral 
whereas larger PAHs are more likely to be charged. Hence, one would expect a mixture of 
neutral and ionized PAHs in the HD 34700 disk. 

Adopting a log-normal size distribution for the PAHs, the HD 34700 disk is in favor 
of smaller PAHs (ao = 3.5 A and a = 0.3) compared to that of the diffuse ISM. This is 
supported by the the prominent 3.3 fin i PAH feature, which is indicative of small, neutral 
PAHs. The 3.43 fin i feature commonly attributed to aliphatic C-H stretch allows us to place 
an upper limit of ~0.22 on the ratio of C atoms in aliphatic sidegroups to that in aromatic 
benzene rings. The 6.89 fin i feature, an aliphatic C-H deformation band, derives the ratio 
of aliphatic C atoms to aromatic C atoms to be ~0.14. 

The total IR luminosity of Z®r « 2.75 L® is derived with a distance of d — 260 pc. This 
leads to Lir/L* ~ 0.13, and the PAH emission accounts for ~ 10% of the total L IR . Our 
best fit model suggests that the dust spatial distribution peaks at ~ 100 AU and decreases 
outward with a power-law index of 7 « 2. Taking the disk geometry into account, the total 
dust mass in the disk is estimated to be Md us t ~ 21.9 M®. Porous grains larger than ~ 30 fim 
attain an equilibrium temperature of T <120 K at r « 25 AU, indicating that ice mantles 
on these grains begin to develop at r & 25 AU. Therefore, r = 25 AU is regarded as the 
“snowline” in the HD 34700 disk. 

Without any prior knowledge of the PAH spatial distribution, the observed PAH emis¬ 
sion yields a total mass of Mpah ~ 4.53 x 10“' (r/ AU) 2 M® if one assumes that all the PAHs 
are at a distance of r from the central star. As free-flying PAHs will easily condense onto 
the ice mantles beyond the snowline, an upper limit of Mpah ~ 2.83 x 10“ 4 M® is derived for 
r = 25 AU. This results in a PAH-to-dust mass ratio of Mpah/ Md us t ~ 1.3 x 10” 5 , indicating 
a deficiency of PAHs in debris disks compared with the diffuse ISM. 

The photodissociation of PAH molecules occurs continuously by absorbing an energetic 
photon, and small PAHs (< 4.1 A) are likely to have been photo-destroyed during the lifetime 
of HD 34700 (r* ~ lOMyr). We calculate the PAH mass replenishment rate to be ~ 4.5 x 
10“ 9 M® yr -1 due to photodissociation. 

The stellar RP dominates the dust and PAH removal in the HD 34700 disk. Due to the 
RP together with Poynting-Robertson drag, ~ 2000 M® of dust has been removed during 
the ~ 10 Myr lifetime of HD 34700. Similarly, ~ 30 M® of PAHs has been removed out of the 
disk. Collisions and/or sublimations of planetesimals, asteroids, and comets are the likely 



origin of the dust and PAH replenishment. 
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